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A problem of the definition of the heat transported in thermomagnetic phenomena has been well 
realized in the late sixties, but not solved up to date. Ignoring this problem, numerous recent 
theories grossly overestimate the thermomagnetic coefficients in strongly interacting systems. Here 
we develop a gauge-invariant microscopic approach, which shows that the heat transfer should 
include the energy of the interaction between electrons and a magnetic field. We also demonstrate 
that the surface currents induced by the magnetic field transfer the charge in the Nernst effect, 
but do not transfer the heat in the Ettingshausen effect. Only with these two modifications of 
the theory, the physically measurable thermomagnetic coefficients satisfy the Onsager relation. We 
critically revised the Gaussian fluctuation model above the superconducting transition and show 
that the gauge invariance uniquely relates thermomagnetic phenomena in the Fermi liquid with the 
particle-hole asymmetry. 



PACS numbers; PACS numbers: 71.10.-w 

I. INTRODUCTION 

It is well known that in the electric field, E = — V0, 
and magnetic field H, the energy of a system of charged 
particles with coordinates and momentunis isi 



^ 2m ^ 



^^(r.xp„).H. (1) 



Considering the heat transfer, it is important to realize 
which terms in the above equation should be associated 
with the thermal energy. The uniform electric field accel- 
erates all particles in the same way, while the effect of the 
magnetic field depends on the particle state. Therefore, 
it is reasonable to assume that the electric term (the sec- 
ond term) contributes to the potential energy, while the 
magnetic term (the third term) provides contribution to 
the thermal energy. Of course, the above consideration 
cannot be considered as a proof. Moreover, there is no 
regular methods for identification of thermal energy. Cor- 
rectness of our assumption will be verified through the 
Onsager relation between the Ettingshausen and Nernst 
coefficients, which describe the heat and electric currents 
initiated by the electric field and temperature gradient, 
correspondingly. We will show that the magnetic term 
overlooked in previous works plays a crucial role and re- 
stores the gauge invariance. While in accordance with 
van Leeuwen's theorem^ the magnetic contribution to the 
heat transfer does not allow a classical interpretation, it 
is critically important for consistent quantum description 
of thermomagnetic transport. 

The definition of the heat current and the Onsager re- 
lation in the magnetic field are long standing problems. 



which were attracted significant attention in the late six- 
ties in the context of thermomagnetic vortex transport 
in the type-II superconductors^''^. Maki^ has noted that 
the phenomenological^ and microscopio^ descriptions of 
thermomagnetic phenomena violate the third law of ther- 
modynamic and the Onsager relation. All attempts to 
resolve this enigma were based on thermodynamic treat- 
ment. Various corrections suggested to the heat current 
are expressed in terms of the equilibrium magnetization 
currents^i^. As we will show, the magnetization cur- 
rents do not transfer the heat and the thermodynamic 
approach cannot solve this problem. 

Recently, much efforts have been dedicated to 
understanding of thermomagnetic effects in high- Tc 
superconductors®'^'*. Variety of phenomenological con- 
cepts were suggested to explain anomalously large Nernst 
signal above T^^ . In many paper a^°i^^'^^i^^ the Nernst ef- 
fect is attributed to the motion of vortices created by 
phase fluctuations. Microscopic modelsi^ii^ are based on 
the Gaussian-Fluctuation Theory (GFT), which includes 
both module and phase fluctuations, but is applicable 
only near Tc^^iiiii^ii^ (for a review see^). However, 
the theories related to vortices are based on the results 
of Refsi^ and^, which inconsistency with basic principles 
was noted vcS^. GFT and other microscopic models also 
ignore corrections to the heat current. 

While the large Nernst coefficient was obtained in 
all above models, there is still not much understand- 
ing of how superconducting fluctuations could lead to 
this effect. In the model of noninteracting electrons, 
which describes well ordinary metals, the Nernst coef- 
ficient is small, because it is proportional to the square 
of the particle-hole asymmetry (PHA)-^^. In this case, the 
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Nernst coefHcient combines PHAs of the thermoelectric, 
r], and Hall, a^y, coefhcients: iV„ ~ rj ■ a^y The inter- 
electron interaction cannot change PHA of rj and (Jx^^- 
According to current points of view^ i-'^'^'-'^^i-'^^i-'^'^i-'^^'-'^^ , the 
interaction in the Cooper channel leads to the Nernst co- 
efficient in the zeroth order in PHA. This effect is (e f /TY 
times larger than the corresponding correction to rj ■ a^y 
(ef - is the Fermi energy). However, experiments with or- 
dinary superconductors, such as Nb, Al, Sn, do not sup- 
port this prediction. This contradiction is a consequence 
of a number of misconceptions at microscopic level. 

The microscopic formalism based on the Kubo method 
allows one to calculate the Ettingshausen coefficient T, 
which describes the heat flux induced by crossed electric 
and magnetic fields, j'* = — T[E x H]. The tempera- 
ture gradient cannot be directly introduced in the La- 
grange formalism of the Kubo approach. Therefore, the 
Nernst coefficient, which relates the electric current with 
the temperature gradient, j*^ = N[VT x H] , is found from 
the Onsager relation, 

N = T/T. (2) 

The coefficient T is calculated for the infinite sample, 
while the Onsager relation is only valid for a finite 
sample^^, where surface currents should be taken into 
account. 

If magnetization depends on temperature, the tem- 
perature gradient generates surface magnetization cur- 
rents, which provide important contribution to the 
charge transfer in the Nernst effect^. Several recent 
worksiSiiSiiSi state that surface magnetization currents 
also play an important role in the heat transfer and the 
Ettingshausen coefRcient calculated for the infinite sam- 
ple, STinf, should be corrected due to the surface heat 
currents: ST = JT„/ - c/x. According toi^ii^iiS, GET 
leads to the bulk coefficient, which is of the same order as 
the correction due to surface currents: STinf = 3/2 c/i, 
where /i is the fluctuation magnetic susceptibility^, 

e^T r 2oi|^ for 2D, 
^ ~ "6^ \ for 3D, 

77 = (T - Tj/Te and a - (T - T,){i{T)f, i{T) is the 
coherence length. 

In the current paper we reconsider basics of the micro- 
scopic description of the thcrmomagnetic effects. This 
work will address the following questions: (1) Is the heat 
current operator modified in the presence of a magnetic 
field? (2) Do surface magnetization currents transfer the 
heat? (3) How do thcrmomagnetic coefficients satisfy 
the Onsager relation? (4) Einally, can the interelectron 
interaction change the PHA of the Nernst and Etting- 
shausen coefficients? While our results are quite general, 
we specify them for GET, which will be critically revised. 
In particular, we show that ST ~ TSN ^ {Td tpY [i- 



II. ETTINGSHAUSEN COEFFICIENT: KUBO 
METHOD 

In the gauge H = i[kjj x Ajj] and E = —i]s.E4>, contri- 
butions of static electric and magnetic fields to the energy 
of a charged particle are easy separated, 

^ ^ (p + cAh/c)^ + e0 = + -(v • Ah) + + /io,(4) 
Zm c 

where = p^/2m — /^O: A'o is the chemical potential. Let 
us first show the importance of the second (magnetic) 
term for noninteracting electrons^^. Thermal energy is 
counted from the electro-chemical potential, so only the 
first and the second terms contribute to the heat current. 
The corresponding diagrams for the Ettingshausen coeffi- 
cient are shown in Eig. 1. The diagram (a) has a form of 
the Hall diagram^^, where the electric current operator, 
ev, is replaced by the heat current, ^pV . The Hall effect 
is proportional to PHA, i.e. the corresponding integrant 
is an odd function of and a nonzero Hall coefficient 
is obtained after expansion of all parameters in ^p/eF 
near the Eermi surface. With the heat current operator 
^pV, the integrant becomes an even function and gives 
a nonzero Ettingshausen coefficient without expansion. 
However, this large contribution is canceled by the dia- 
gram (b). The well-known expression for T in a system 
of noninterating electrons is obtained in the second order 
in PHA, when both integrants (a) and (b) are expanded 
in (^p/ef)^ (for details see an Appendix in2^). 

Now we consider electrons interacting in the Cooper 
channel. The energy current operator in this system^! is 
easily expanded to include external fields, 

J' = ^ v^pfl+flp + ^ — (vAi/)a+ap + ^ ve^a+Op 
p p p 

+ ^ a+flp - A/2 XI + ap+p'-p"«p"«p'«P 
p p,p',p" 

+ X (v + v')C/i„ip(Ri) exp[z(p + p')i?i] a+flp/ , (5) 

p,p',R, 

here a+ and Op are the electron creation and annihilation 
operators, A is the interaction constant in the Cooper 
channel, and Uimp is the impurity potential. The first 
and the forth terms in the Eq. IC2I describes the energy 
flux of noninteracting electrons without flelds. The sec- 
ond and third terms are due to electron interaction with 
magnetic and electric flelds. Two last terms are due to 
the electron-electron and electron-impurity interactions. 
These two terms generate diagram blocks (Aslamazov- 
Larkin blocks) proportional to PHA, but contributions 
of all blocks with PHA cancel each other^— . 

The thermal energy is defined as the electron energy 
counted from the electro-chemical potential e4> + fj, and, 
therefore, the third and the forth terms in Eq. IC2l do not 
contribute to the heat current. Thus, calculating the heat 
current we should take into account only two heat current 
vertices 7^ and 72 corresponding to the first (kinetic) and 
second (magnetic) terms in Eq. IC2I 



3 



Two leading terms in the heat current operator gen- 
erate two diagrams, which describe T in the Aslamazov- 
Larkin (AL) approximation (see Fig. 2). The wavy hnes 
correspond to the fluctuation propagato r^°i^^ , 



-R,A 



(6) 



2 CjLUd 



- ag^ ± — + 7^ ) ,(7) 



where P{q, lo) is the polarization operator, v is the elec- 
tron density of states, lud is the Debye frequency, and 
C-y is the Euler constant. The last term in Eq. [7] is 
proportional to PHA^ii^. 

The AL blocks B^'^'^ presented in Fig. 2 are built 
from electron Green functions and vertices j'^-'^'^ (see 
Tab. 1). The left blocks B^* and B^ in Figs. 2. a and 
2.b are blocks with heat current vertices 7^' (kinetic) 
and 72 (magnetic). The right block B*^ in both dia- 
grams includes the electric current vertex 7*^ = ev ■ e^;, 
Be — E/E. Block B^ includes the magnetic vertex 
7^ = {e / c)v ■ Ah ■ Results of calculation are summarized 
in Table L AL blocks are obtained by insertion vertices 
7 into the polarization operator and can be expressed 
through VqP^(q,0). Blocks B% B^, and B^ are well 
known^'^. The block B'j* describing the heat current in 
the absence of magnetic field has been calculated in22 
(see alsoiSii^). Here we introduce B2, which is based on 
the electron vertex 72 and describes the magnetic correc- 
tion to the heat current. 

The first AL diagram (Fig. 2. a) was investigated in 
Refji^ and its contribution is ^T^^^ = 3/2 c^. The same 
result has been obtained in the time dependent Ginzburg- 
Landau formalism (TDGL)i^ii^i^. 

Contribution of the second AL diagram (Fig. 2.b) is 



"(-(2) 7J _ CS, 



dq duj B^B'^ 



(27r)" 2tt 2n 



{L%L^+L'lL'l), (8) 



where = coth(tj/2r)(L^ - i^), L± is used for 
L(q ± k/2, u} ± Vl/2), and n is the system dimensionality 
with respect to the coherence length ^(T) . Expanding the 
integrant to the linear order in and k and calculating 
the integrals over oj and q, we find that the contribu- 
tion of the second diagram, t|^^^ cancels completely the 
contribution of the first one. 

Thus, without PHA the Ettingshausen effect is absent, 
(5T = 0. To get nonzero result, we should expand the 
fluctuation propagator (Eq. [B]) up to the second order 
in PHA. Expanding the polarization operator (Eq. [7|) to 
the second order in 70; we get 



T Attc 

where 7 = — In — - — — Thus, thermomag- 

2eF omeF ttTc 
netic coefficients in the fluctuation region are propor- 
tional to {T/tFY- 



r 2aM 



for 2D, 
for 3D; 



(9) 



1 d\\iv 2C^ujd 



Summarizing results of this section, we would like to 
note that the total heat current operator of the fluctuat- 
ing pairs in the magnetic field. 



B^ 



Bj + B^ 



[q+(2e/c)Aff], 



(10) 



may be considered as the gauge-invariant extension of 
the operator B^* without H. This is a key point, because 
the further calculations of diagrams in Fig. 2 are similar 
to that for noninteracting electrons (Fig. 1): the kinetic 
and magnetic terms in B'' generate two diagrams, which 
cancel each other in the zeroth order in PHA. 

In the above calculations the interaction with the mag- 
netic field has been included in the heat current. Note 
that Eq. [TO] can also be derived in another, more formal 
approach, where the magnetic field is initially included 
in electron states. Without magnetic field, the thermo- 
electric coefficient is described by the AL diagram with 
the heat and electric current operators, Bf' and B'i^iSl. 
In the magnetic field, the momentum of the Cooper pair 
is given by q + 2eA/f /c, and the polarization operator 
has a form: P(q+2eAij/c, uj)'^'^. Calculating the ther- 
momagnetic response, all blocks of the diagram should 
be expanded in Ah- Expanding polarization operators 
in the fluctuation propagators, we obtain the flrst dia- 
gram for the thermomagnetic coef5cien t^^i^° . Expanding 
the heat current block B^* = w VqP^(q -I- 2eA///c), we 
immediately obtain the block Bj, which forms the sec- 
ond diagram. This magnetic term has been lost in all 
previous work o-'^^'^^i^^'^^ . Using Eq. 10 as the heat cur- 
rent operator for fluctuating pairs, our results can be also 
obtained in TDGL formalism. 



III. NERNST COEFFICIENT: QUANTUM 
TRANSPORT EQUATION 

To investigate a response of the electron system to VT, 
one can use the quantum transport equation. Previously 
we adapted this method to calculations of the themoelec- 
tric and Hall coefficients in GFT^^'^^. In this approach 
the electric current is given by 
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dq duj 



(27r)" (27r) 



qlm <5L^(q,w), 



(11) 



where SL'-' (q, oj) is the nonequilibrium correction to 
the fluctuation propagator. In the equilibrium, — 
IjRpC i^A^ where the Keldysh component of the polar- 
ization operator = iirv/A at T - Tc < Tc. 

Calculation of (5L'-^(q, w) for the Nernst coefficient is 
analogous to its calculation for the Hall effect^. The 
nonequilibrium effects are taken into account by the 
VT and H-Poisson brackets between polarization oper- 
ators: {Pi, Pair = yTidrPi ■ 9,P2 - dTP2 ■ dqP{) and 
{Pi,P2}h = (e/c)H- [a„Pi X 9,P2], where drP = dP/dT 
and dqP = dP/dc^^^L21. Therefore, the only differ- 
ence in calculations of the Hall and Nernst coefficients is 
that the derivatives dP/doj in Hall coefficient (E-Poisson 
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bracket22.) should be replaced by the derivatives dP/dT 
in the Nernst coefficient. 

To get the vector product H x VT, the H-bracket 
should include the same polarization operator as the 
VT-bracket. Thus, in the first order in H x VT the 
nonequilibrium fluctuation propagator (5L'^(q, cj) is de- 
scribed by the diagrams shown in Fig. 3. Finally, calcu- 
lating 5L'~^{(\,uj) and the corresponding Nernst current 
in the interior of the sample (Eq. Ilip . we find 

where 5N is equal to the term, which was calculated in 
the previous section from the Onsager relation (see Eq. 

Thus, in the infinite sample or in the interior of the 
finite sample, the Nernst coefficient consists of two terms. 
The second term in Eq. [T^] has the zeroth order in PHA 
and violates the Onsager relation. As it will be shown in 
the next section, in the finite sample this term is canceled 
by the contribution of the magnetization currents. 



IV. ONSAGER RELATION IN MAGNETIC 
FIELD 

The above results, Eqs. [9] and [121 have been cal- 
culated for the infinite sample. Referring to^, recent 
worksiSiiSii^ state that for a finite sample both coeffi- 
cients should be corrected due to charge and heat trans- 
fer by surface magnetization currents. Here we show that 
the magnetization currents contribute only to the charge 
transfer, and the results o^ have been misinterpreted. 

The electric magnetization current y^^g in the poten- 
tial relief 0(r) transfers the energy fiux j^^^ = cjy'f^^g 
(Eq. 37 in?3). Using = cfiPAn, we get 



is associated with the bulk current in the finite sample, 

N^nf = itnr/NT X H] . Usiug Eq. M WC get 



Jmag 



= cfi[li X E]. 



(13) 



This term is erroneously attributed to the heat 
fiux^§.ii2ii^i2i. As we discussed, the electric potential cj) 
and the corresponding vertex 7"^ — evcj) do not contribute 
to the heat current, because the thermal energy should 
be counted from the electro-chemical potential. 

In the interior of the sample the electric current con- 
sists of the transport and magnetization components, 
iinr = itr + Jman- "^^^ magnetization component is^ 



c||(VTxH) 



(14) 



The magnetization currents are divergence-free. The 
total magnetization current through the sample cross- 
section must be zero, i.e. the bulk magnetization cur- 
rents are canceled by the surface currents. Therefore, the 
Nernst coefficient measured in the finite sample is deter- 
mined by the transport currents: N = jj,./[Vr x H] -"^. 
The Nernst coefficient in the infinite sample (Eq. [1/ 



SN = 



•itnag 



[VT X H] 



= mnJ-C^. (15) 



Taking into account Eq. [3l we see that the second term in 
the last equation, cdfi/dT, cancels completely the second 
term in SNinf (Eq. [T2|) . The rest is equal to SN, which 
satisfies the Onsager relation SN = ST/T. 

Note that, if in contradiction to our results, the surface 
magnetization currents provide the heat transfer, this ef- 
fect could be found in transport measurements. In the 
Gorbino disc geometry with the magnetic field perpen- 
dicular to the disc and the circular inductive electric field 
in the plane (Fig. 4), the heat current in the radial di- 
rection does not contain the surface components, which 
were predicted for the standard parallelepiped geometry 
^j ^i8,i9,24 ^ According to our results, both experiments 
will give the same results. The surface electric currents 
generated by VT are very significant (see Eq. [Tl)) . How- 
ever, they cannot be experimentally separated from the 
interior currents, because contrary to the circular elec- 
tric field, the circular temperature gradient does not ex- 
ist. This difference between E and VT is refiected in the 
asymmetry of the Nernst and Ettingshausen coefficients 
calculated for the infinite sample. 



V. PARTICLE-HOLE ASYMMETRY IN 
THERMOMAGNETIC EFFECTS 

Now we show that in the general case, the interelec- 
tron interaction cannot change PHA requirements for N 
and T, i.e. the thermomagnetic coefficients are always 
proportional to the square of PHA. 

Assuming that electron scattering from impurities is 
the main mechanism of the momentum relaxation, it is 
easy to see^^i^li^ that the magnetic field and tempera- 
ture gradient enter into the transport equation formal- 
ism through the distribution functions of noninteracting 
electrons and the Poison brackets. In fact, the terms 
proportional to VT x H can appear in three different 
ways^^: (a) through the Nernst nonequilibrium distribu- 
tion function of noninteracting electrons, (er^/cm) v ■ 
[VT X H] (dS/dT), (b) through the H-Poisson bracket 
that involves the nonequilibrium distribution function 
under the temperature gradient, — eT(v-E)(9S'/i9T), and, 
finally, (c) due to double, VT and H, Poisson brackets. It 
is evident that a and &-type terms have already includes 
the Hall PHA, which is proportional to {dv/dp) = 1/ni. 
The c-type terms in the form of the double Poisson brack- 
ets describe the AL process, which has been investigated 
above. As we have seen, the AL diagram gives the con- 
tribution in the zeroth order in PHA, however, this con- 
tribution is canceled by the contribution of the surface 
magnetization currents. Thus, we conclude that the in- 
terelectron interaction can provide many-body thermo- 
magnetic effects only in the second order in PHA. 
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VI. CONCLUSIONS 

We have shown that the magnetic term in the Hamil- 
tonian of charged particles (Eq. [T]) should be associated 
with the thermal energy. The corresponding term in 
the heat current operator (Eq. 2]) restores the gauge in- 
variance and gives important contribution to the Etting- 
shausen coefficient. We also found that the surface mag- 
netization currents do not contribute to the heat current, 
but provide substantial contribution to the charge trans- 
fer in the Nernst effect (Eq. [T5)) . Our gauge- invariant 
scheme gives the thermomagnetic coefficients that sat- 
isfy to the Onsager relation (see Eqs. [15] and [T^ . In 
the general case of the Fermi liquid with particle-hole 
excitations, we conclude that the measured thermomag- 
netic coefBcients are always proportional to the square of 
PHA. Any interaction by itself, i.e. without changing the 
electron band structure or character of elementary exci- 
tations, cannot provide large thermomagnetic effects. 

The developed approach has been applied to effects of 
superconducting fluctuations. We show that the gauge 
invariant form of the heat current operator of fluctuat- 
ing pairs is i/aa;(q -I- 2eAi//c). The second (magnetic) 
term missed in previous publications plays an important 
role: as in the case of noninteracting electrons, the cor- 
responding diagram cancels completely the large, zeroth- 
order in PHA term in the Ettingshausen coefficient. We 
also show that the Nernst coefficient in the infinite su- 
perconducting sample consists of an anomalously large, 
the zeroth order in PHA term (Eg. [T^ . However, in the 
finite sample, this term is cancelled by the the surface 
magnetization currents (Eq. llSp and thermoelectric co- 
efficients satisfy the Onsager relation. Our results for 5N 
are different by a huge factor of (ep/T)'^ from the pre- 
vious works^^JiiiSiii^, which claim that the attractive in- 
teraction in the Cooper channel provide thermomagnetic 
transport without PHA at all. We can also easily rebuff 
this claim, if we just change a sign of the interaction con- 
stant and consider the repulsive interaction in the Cooper 
channel. As known, in ordinary metallic films this inter- 
action results in corrections to conductivity, which are a 
factor of (epr)^^ smaller than the conductivity of non- 
interating electrons'^^ . If the statement of Refs i^^i^^'^^i^^ 
is correct, the thermomagnetic effects in ordinary metals 
would be {eFT)-\eF/Tf - {eF/T)/{TT) larger than 
predicted for noninteracting electrons. Certainly this 
huge effect is not known. 

Returning to the problem of high- Tc superconduc- 
tors, we should note that most of the models, including 
fluctuation exchange^^, 'cheap' vortices^, and preformed 



Using GET as an example, we have shown that the 
gauge-invariant form of the heat current is critical for 
description of the Ettingshausen effect and that the sur- 
face currents are important for the Nernst effect. All 
other models of thermomagnetic transport including vor- 
tex models should be reconsidered in accordance with the 
formalism developed above. 

We would like to acknowledge useful discussions with 
I. Aleiner, A. Larkin, D. Livanov, A. Varlamov, and I. 
Ussishkin. 



APPENDIX A: METHOD KUBO FOR 
NONINTERACTING ELECTRONS 

Appendix A 

In this appendix we present detailed calculations of 
the Ettingshausen coefficient for noninteracting electrons 
using the Kubo method. 

According to Eq. 4, the heat current vertex for nonin- 
teracting electrons in the magnetic field is given by 



l"" = li + 72 = ?pv + (e/c)(v • A)hv. 



(Al) 



Two vortices 7^ and 72 create two diagrams shown in 
Fig. 1. Solid lines in diagrams represent the electron 
Green functions. 



(A2) 



To get the Ettingshausen coefficient proportional to A(k- 
E), one should expand the Green function G{p + k) in 
powers of (k • v). Then the contribution of the first dia- 
gram is given by 



cH 



dp de dSo 



Cpv(vA)(vE)(vk)/i,(A3) 



where Sq — — tanh(e/2r), n is the dimensionality of the 
system, and the combination of the Green functions is 

h = 2iGfG^lm{Gff = 2iTHm{G^f. 

The contribution of the second diagram is 



(A4) 



To = ^ 



dp de dS( 



v(v A)(vE)(vk)/2, (A5) 



cH J (27r)" 27r de 

where the combination of the Green functions I3 is 

I2 = 2iGpm{G^f = 2iT^lmGp (A6) 

Integration over angles of the electron momentum in Eqs. 
(A. 3) and (A. 5) gives 
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ion excnange , cneap vorrices , ana preiormea /• yi 

have treated the state with large Nernst effect as / d^lp v(v • A)(v k)(v • E) = —— — A(k • E). (A7) 



pair; 

natural extensions of GET. In the light of our results, the 
concept of this specific state which continuously evolved 
from GFT due to fluctuations and additional interactions 
is far from obvious. In this sense, the current work rather 
supports an idea of the Nernst state as " a thermodynamic 
phase in itself with characteristic properties speciflc to it 
as a fixed point, distinct from those of a normal metal"—. 



Then the total contribution may be presented as 
"2 de 95-0(6^ 



T = 



n(n + 2)cHj 2tt'^^^ de 



A(k • E) 

X v''T'v[£,plm{G^f + lmG^]. (A8) 
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Without taking into account PHA the total contribution 
of two diagrams goes to zero after integration over ^p, 
because 



0. 



(A9) 



Nonzero contribution arises from terms proportional to 
e^, thus we should expand all electron parameters near 
the Fermi surface. For example for 3D conductor, 



4 
V V 



2 ep 



, (AlO) 
(All) 



Taking into account terms proportional to the square of 
PHA, e.g. ^^/e| or ^e/ej,, we get 



,5 

4 £2 



2 m ep 



(A12) 



Substituting this result into Eq. (A. 8), and performing 
integration over e, we get the well-known result for the 
Ettingshausen coefficient of noninteracting electrons. 



31? 



f^To — 

b ep H 



(A13) 



where — eH/mc is the cyclotron frequency and Uxx 
is the Drude conductivity. For 2D-conductors, the cor- 
responding relation between and two-dimensional 
conductivity has an additional numeric factor of 2. 

Thus, calculating the Ettingshausen coefficient of non- 
interacting electrons, we demonstrated that the magnetic 
field should be taken into account in the heat current ver- 
tex. The diagram with this magnetic vertex in the heat 
current cancels the basic diagrams in the zeroth order in 
PHA. The nonzero Ettingshausen coefficient arises only 
in the second order in PHA. As it has been shown in the 
main text, the above conclusions are also relevant to any 
many-body corrections to thermomagnetic coefficients. 



The block B*^ with the electric current vertex, 7*^ = 
ev • e^;, may be presented as^ 

B^(q) = 2e VqP^(q,0) ■ es = 2ei/aq-es. (B3) 

The block with the vertex 7-^ = (e/c)v • Ah is given 
by20 



2e 
c 



(B4) 



The block Bj' with the kinetic heat current vertex, 
7? = • (eyx = ||A), is given bj^ (see 

alsoi^^i^O) 

Bi(q,w) = Lo VqP^(q,0) • ey. = uva q • ey.. (B5) 

Next, we calculate the block B2 with the magnetic heat 
current vertex 72 = (v • A//)(v • e^h). The integral over 
angles of the electron momentum involves only the vertex 
72 , because the heat current is in the direction of Ah- 
To obtain an imaginary part in Eq. IBll the integral 



dS, (Gp ) Gq-p 



27ri 



(2e LU q - V i/r)' 



r, (B6) 



should be expanded in lu (in calculations of Be it is ex- 
panded in q ■ v) . Finally, we get 

B^(q,c^) = 2toVlP"iq,0)AH = 2(e/c)wi^aAff . (B7) 



APPENDIX C: GAUGE E = ifiAjs/c 

Here we will show how our results can be obtained in 
the gauge, where E = iilAp/c and H — i[k x Ah] (this 
was a question of one of our referees) . In the electric and 
magnetic fields the kinetic energy has a form K ~ {p + 
eA/c)2/2rn, and the part of the Hamiltonian describing 
the interaction with external fields is given by 



APPENDIX B: AL BLOCKS 

For an arbitrary electron momentum relaxation time r, 
the AL blocks B'^''''^ built from electron Green functions 
Q^-i^) with electron vertices 7 (7*^, 7^, Ji, and 72 ) are 
given by2£ 

Tie,h,H J j dp de f,fiH n I \ iGp)'^G^_ 
B,' =Im / T^-^— 7,- So{e) — f (Bl) 



C 



(27r)" 27r 



(B2) 



TTiyr J (27r)3 ^ 
where the electron Green functions are given by Eq. A. 2. 



Calculating a response to E x H = (ri/c)[A/f (k • A^;) — 
\s.{Ae-Ah], it is convenient to use the gauge conditions k- 
A = and Ah ■ Ap = 0^^ . In this gauge the second term 
in Eq. ICll can be neglected. Including the interaction 
with the magnetic field we get the heat current operator. 



J' = «P«p+E v^^'^^^ ^^^^ 



As it is expected, the term in the heat current describing 
the interaction with the magnetic field is independent on 
the presentation of the electric field (see Eqs. 5 and Al). 
Therefore, all further calculations of thermomagnetic co- 
efficient are the same as in the gauge E = —V(j>. 
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TABLE I: AL blocks (operators for fluctuating pairs) B based on electron operators 7. 



1" 


= e V • cb 


7^ = (e/c) V • A 


7i 


= ^ V • e^i. 


h 
72 


= (V • Aif)(v • Bjh) 


B*^ 


= 2e VqP'^Cq, 0) • be 


B" = (2e/c) VqP''(q,0)-Ajy 


B? 


= uj VqP^(q, 0) • e^h 


B^ 


= 2uj V^P'*(q,0) Ah 




= 2eva q • es 


= (2e/c)z/a q • A/f 




= wi^a q • Gjh 




= (2e/c)wz/a Ajf • e^h 



Av 

^pv/^ ^ev (Av)vC 7ev 



Fig. 1. Diagrams for Ettingshausen coefficient of noninteracting electrons. 
The straight lines stand for the electron Green functions. 




Fig. 2. Fluctuation AL diagrams describing the heat current in crossed electric and 

magnetic fields. Wavy lines stand for the fluctuation propagators and straight 
lines stand for the electron Green functions, which form the AL blocks. 

vr 




Fig. 3. Nonequilibrium fluctuation propagator SL due to H and VT - Poisson brackets. 
I 



I 




Fig. 4. The interior heat currents in the Ettingshausen effect can be measured in 
the Gorbino disc geometry with the circular inductive electric field. 
The interior electric currents in the Nemst effect cannot be measured, because 
the circular temperature gradient does not exist. 



